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Abstract 

A dynamic response to a magnetic field in a long disordered cylinder is 
considered. We show that, although at high frequencies conduction is classical 
in all directions, the low frequency behavior corresponds to localization in the 
longitudinal direction and to a diamagnetic dynamic persistent current in the 
transversal one. The current density does not vanish even in the limit of the 
infinitely long cylinder. Being of a purely dynamic origin the current can be 
destroyed by inelastic scattering but at low temperatures the decay time can 
be very large. 
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The recent measurement of the magnetization of an array of 10 7 isolated Cu-rings 
attracted a lot of attention to the problem of persistent currents in disordered systems. The 
circumference of these rings was much larger than the elastic mean free path and therefore 
at short distances electron motion was diffusive. Considerable effort has been taken in order 
to build a proper theory of the observed effect in the diffusive regime. Calculation of the 
persistent current by averaging over disorder at fixed Fermi energy gives in this regime an 
exponentially small amplitude. In order to obtain a non vanishing effect one should either 
take into consideration electron-electron interactions |§ or calculate with a fixed number of 
particles In both cases it is essential that all the sizes of the sample are small. It 

is clear that for a macroscopic object the difference between results for the canonical and 
grand canonical ensembles cannot be important and therefore one should not hope to have 
a thermodynamically stable persistent current in such a case. 

However, does it mean that there is no possibility to have a long living macroscopic 
current in a macroscopic piece of a disordered metal? Of course, if such a current did not 
correspond to the minimum of the free energy different types of inelastic scattering would 
kill the effect. Nevertheless, the corresponding decay times would be very large and strongly 
dependent on temperature which makes a discussion of such a possibility quite interesting. 

In order to understand better what type of effect we want to consider below let us 
imagine first a very long clean metallic cylinder with the length L z ,the circumference L and 
the thickness of the walls d. Applying a magnetic field H along the cylinder one induces 
currents. In the ballistic regime the cylinder is an ideal diamagnet and the induced total 
current I d is equal to 

I d = dL z LK d H/4nc, K d = -ne 2 /m. (1) 

In Eq.(|l|) m is the electron mass, e is the electron charge, n is the electron density, c is the 
light velocity. The current Id can exist at very low temperatures for a very long time but it 
does not correspond to the free energy minimum and must finally decay. Experimentally, the 
current can be obtained first cooling down the cylinder and then applying the magnetic field. 
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The inverse procedure would lead to a thermodynamic current which can be much smaller. 
The averaged thermodynamic persistent current I t h for clean cylinders was calculated in ref. 
H and can be estimated as 

hh ~ \h\ /(poL'f 2 , (2) 

where po is the Fermi momentum. Eq.(|2]) shows that the thermodynamic persistent current 
can be much smaller than the dynamic one which, of course, is not unexpected. The dif- 
ference between the dynamic and thermodynamic currents is due to existence of degenerate 
states. If one represents energy levels by functions of the magnetic field the degeneracies 
correspond to level crossing. In the diffusive regime Pq 1 <C I <C L, where I is the mean free 
path, the thermodynamic current is proportional to exp(— L/l) || and is negligibly small. 
In this Letter we calculate the corresponding dynamic current and show that it can be quite 
noticeable. For simplicity we consider a standard model of non interacting electrons moving 
in a random potential. Calculating the dynamic current we use standard formulae for linear 
response. We assume that the external magnetic field H has the form 

H = H + H u cos{ut), (3) 

where H is the static component which can be arbitrary. The second term in Eq.(||) describes 
the oscillating part of the magnetic field which induces an electric field directed along the 
circumference and, correspondingly, an oscillating current. We assume that the amplitude 
H u of the oscillating part of the magnetic field is small and calculate the induced current 
using the conventional linear response theory. The density of the oscillating current can 
be written as 

j u = c^K^Au, (4) 

where K (u) is the current-current correlation function and A u is the vector potential cor- 
responding to the magnetic field in Eq.(|3|). We use the gauge with the vector potential 
directed along the circumference. Using the retarded Gf fr, r'j and advanced {r, r'^j 
Green functions we can write the function K (u) as follows 
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K (u) = Kt (co) + K 2 (u) , 



(5) 



K 1 (u) = (e 2 /iTc) [ [n(e-^-n(e-u-^](v r Gf(ry)v r >Gt^(r\r))dr'de 1 (6) 



K 2 {u) = (e 2 /iir) fn(e- v r Gf +[jJ (r,r'j tyGf (r',r) -v r Gf (r, r') v r >G^_ w (r',r 



dr de+ 



(e 2 /iirm) J n(e — //) ^ Gf (r, r 



G A (r, r) 



de. 



(7) 



where v r = (1/m) [— zV — (e/c) A] is the component of the velocity operator perpendicular 
to the axis of the cylinder, A is the vector potential corresponding to the static component 
H of the magnetic field, n (e) is the Fermi distribution, and the angular brackets stand for 
averaging over impurities. The terms K\ (u) and K 2 (u>) have completely different origins. 
In the limit uj —>■ the term K 2 (u>) can be obtained from the free energy by differentiating 
at a fixed chemical potential over the vector potential and, as we have mentioned, can be 
much smaller than the first one. In the diffusive regime the average over impurities of a 
function containing a product of only retarded or only advanced Green functions can be 
substituted by the product of the corresponding averaged Green functions (^G 11 '^ . In the 
momentum representation this is given by 



(G^ A ) = {e-e{P)±^/2r)- l 



(8) 



where r is the mean free time. Substituting Eq.@ into Eq.(|7|) and integrating the first term 
in Eq. (|7D over momentum by parts one can see that in the limit u — > it cancels the second 
one. Here, we consider the situation when the magnetic field is concentrated in the cavity 
of the cylinder. Otherwise, we would have to sum over the Landau levels and would obtain 
that the quantity K 2 gives additionally the Landau diamagnetism. 

In the following we concentrate on calculation of the quantity K\ (ou) . Before carrying 
out explicit calculations in the diffusive regime let us demonstrate how one can calculate 
this quantity in the ballistic one. Using the representation 



and taking the limit w->Owe obtain from Eq.@ 

\xmK 1 (u) = -(2e 2 /V) £ \v aP \ 2 5 (e a - y) , (10) 

uj— »U — 

where /i is the chemical potential, is the volume. 

In the clean macroscopic cylinder under consideration the eigenstates are plane waves 
and the matrix elements in Eq. fllOD can easily be calculated. All sums are to be substituted 
by an integral over momenta and we can see that lim^o K\ (u) = Kd . 

The classical limit in the diffusive regime corresponds to high frequencies u. In this limit 
one can substitute the Green functions in Eq. @ by their averages Eq. (§) and then calculate 
the integral over momenta. As the result one has 

Kx{u) = -iw(jQ, (11) 

where <To = e 2 Tn/m is the Drude conductivity. 

The limit of low frequencies requires more sophisticated methods of calculation. In 
the diffusive regime the supersymmetry method developed by one of the authors |7|] is 
by now the only possibility to obtain results analytically. This method has been used 
to study localization in long disordered wires |7|||. As we will see below the problem of 
localization is closely related to HJT0| the problem of the dynamic persistent currents. The 
responses both in the longitudinal and transversal directions can be reduced to calculation 
of functional integrals over supermatrices Q. Calculating a correlation function for a model 
of a wire corresponding to a longitudinal response for the cylinders it was shown in Refs. 
ITJII that the system has dielectric properties in the longitudinal direction. Here, we present 
a calculation of the transversal response describing the persistent currents. In Refs. |p|JT0[ 



a general expression for the function K\ (u) has been presented in terms of integrals over 
supermatrices Q with a free energy functional of the non linear supersymmetric a -model. 
This expression has the following structure 



K 1 {u) = -iu(B{Q)) Q , (12) 

where B (Q) is a sum of products of elements of the supermatrix Q. 
In Eq. (|T2| ) the notation (. . .)g stands for the functional average 

(■■■>0 = / (...)ex V (-F[Q])DQ, (13) 

where F [Q] is the free energy functional of the non-linear supersymmetric cr-model. 

The form of the function R (u) which denoted in Refs. flJTCfl the average {B {Q))q remains 
unchanged for the sytem studied, here. While the dynamic response for the problem of 
persistent currents in mesoscopic rings is calculated using the zero dimensional version of the 
cr-model P,|T0||, for the problem of localization in wires one has to take the one dimensional 



version of the cr-model f^JU- Alhtough for the cylinder in the diffusive regime the one 
dimensional cr-model is valid, let us consider a more general model of a stack of small metallic 
rings. Changing the probability of tunneling from ring to ring we can describe the crossover 
from the case of the homogeneously weakly disordered cylinder to the case of isolated rings. 
For such a model the derivatives in the longitudinal coordinate are substituted by finite 
differences. Then the free energy F [Q] takes the form 

F[Q] = - £ J^STr [Q^] + £ F [Q t ] , (14) 



e ^ 2 



-V ( - A[Qi, r 3 ] ) + 2iuKQi 



(15) 



where S = Ld is the cross section of the cylinder, h is the height of each ring, i and j 
enumerate the rings in the stack, Dq is the classical diffusion coefficient. ( We remind that 
Q 2 = 1). The function F [Q] in Eqs.fll4l,l5D describes electron motion inside the rings, the 



first term in Eq.(|l4]) stands for coupling between the rings. We assume that only nearest 
neighbors interact, so that = J for the neighbors, and = otherwise. The limit J = 
corresponds to the stack of the isolated rings. In the limit J ^> 1 the model on the lattice 
Eq.(|i~4"D can be substituted by the continuous one and we return to the one dimensional 



cr-model. The notation STr stands for the Supertrace, the definition of the matrices T3 and 
A can be found in previous works J^Jfi|. 

To recover the case of the homogeneously weakly disordered cylinder one should relate 
J to Do as 

J = nuSD /8h. (16) 
The main contribution for the transversal current when calculating the integral Eqs.fll^, |13D 



comes from terms containing products of two matrix elements of Q taken at the same point 
( terms R\ (u) and R3 (u) in Refs. PJlOfl) and the substitution of the weakly homogeneous 
cylinder by the stack of the rings does not change the form of these terms. 

Due to the one dimensionality of the model Eq.fll4|) one can use the transfer matrix 
technique. Corresponding partial differential equations have been written in Refs. [|7||| for 
the continuous cr-model. For the model on the lattice, Eq.([L4]), the corresponding recurrence 
equation is an integral. Analogous integral equations were written in Refs. ||TT| , |T2|| for the 



model on the Bethe the lattice and also for one dimensional chains. Repeating the main 
steps of these works we reduce calculation of integrals over Qi for all cites % to one integral 
over Q. The corresponding changes can be done substituting Eq.([T3|) by 

(. . ) Q = J (. . .) vl/ 2 (Q) exp (-F Q (Q)) dQ. (17) 

The function \1> (Q) in Eq.([P7D satisfies the equation 

* (Q) = J exp(2JSTrQQ') exp (-F (q')) * (q') dQ' . (18) 

The solution of the Eq.(|18D depends on the vector potential A entering Fq Eq. (p!5|) . In 
principle, it can be solved, at least numerically, for arbitrary magnetic fields using the pa- 
rameterization proposed recently |13| but we will present results only in the limits of zero 
and high magnetic fields. These limits correspond to the orthogonal and unitary ensembles. 
In both cases one can omit the first term in Eq.flT5D because in the unitary case the super- 
matrix Q commutes with the matrix T3 [0. Due to the symmetry of the free energy the 



function \1/ (Q) depends only on the variables A and Ai,2 parameterizing the supermatrix Q 
0. Therefore, in the integral over Q in Eq.(^) one can integrate first over all other variables 
reducing thus the integral to an integral over A and Ai^- At high frequencies deviations of 
the supermatrix Q from A are small and one comes to Eq. ([TT]) . In this limit the conduction 
is classical in all directions. 

Calculations for arbitrary frequency are most simple for the unitary ensemble. Corre- 
sponding calculations are not very different from those performed in Refs. [||[10[ and we 
obtain 

K™ tt (to) = -iuao jl + i J 1 ^°°exp [{mu/A - 5) (A x - A)] ^ 2 (A, Ai) d\d\^ , (19) 

where A = (uSh)~ 1 is the mean level spacing in one ring in the stack and 5 — » 0. 

Non zero persistent currents are possible if at low frequencies u the second term in 
the brackets in Eq. (|T9|) is proportional to Xjio. In the most interesting limit uo — > all 
calculations when solving Eq.(|18D become more simple because the main contribution in all 
integrals over Aicomes from Ai ~ Then the function $ (Q)depends only on one variable 
z = 2u;Aiand we obtain for the response K\ (u) 

Kr U (0) = -a A eff (J) /tt, (20) 

where 

/•OD 

A eff (J) = A exp {-z) V 2 j (z) dz } (21) 
Jo 

is a non trivial function of the coupling J between the rings. This function is known numer- 
ically for arbitrary J. [0 In the limit J = the function ^ = 1 and 

A eff (0) = A. (22) 

In the opposite limit J ^> 1 the function \1/ (z) is the solution of the differential equation 

zd 2 ^/dz 2 - = 0. (23) 

Solving Eq.fl2~3D and calculating the integral Eq. (pl|) we find 
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A e// (J) = A/ (96 J) . (24) 

As we have mentioned above one has localization in the longitudinal direction. In the unitary 
ensemble the localization length L c was calculated in Ref. @ and can be related with the 
help of Eq.(|16"D to the coupling J as L c = 16 Jh. In this limit the function A e jy (J) can be 
rewritten as 

A e// (J > 1) = (6ZASL,)- 1 . (25) 



We see from Eq. fl25|) that in the limit of large J the response Kf nit (0) Eq. (p0|) looks as 
if the cylinder consisted of rings with the height of the order of L c . Compairing Eqs. ([□] 
and |2(|) we see that the characteristic frequency of the crossover from the quantum to the 
classical regime is of the order A e jy. In the model of the weakly homogeneously disordered 
cylinder we can see with the help of Eq. ([16]) that the response Kf mt (0) does not depend on 
the disorder. In this limit it is very small. Decreasing the coupling J makes the localization 
length shorter which leads to a larger value of the response. 

Analogous calculations can be carried out for the orthogonal ensemble corresponding to 
the zero static component H of the magnetic field Eq. (|3|). In this case we obtain as in Refs. 
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K° rm (0) = 0. (26) 

Changing the static component H of the magnetic field we can have a crossover be- 
tween Eqs. (|26|) and (|20|). The characteristic flux C of this crossover is of the order of 
0o(A e ///-E c ) 1/ ' 2 , where E c = ir 2 D /L 2 is the Thouless energy, O is the flux quantum. The 
whole dependence of the response on the flux is periodic with the period 0o/2. 



Adding magnetic or spin-orbit impurities changes Eqs. (|2(]. |26|). The system with the 
magnetic impurities corresponds to the unitary ensemble. One can use as before Eq. ( |20| ) 
provided A e // (J) is substituted by A e j/ (2 J) /2. If the magnetic impurities are absent spin- 
orbit ones lead to a different function A e // (J). However, this difference is only numerical 
and does not change the sign of the response which is in all cases diamagnetic. 
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In conclusion, we showed that a magnetic field applied to a disordered cylinder parallel 
to the axis induces a macroscopic diamagnetic current. In the absence of inelastic scattering 
which can correspond to low temperatures this current can live for a very long time. The 
current density remains finite even in the limit of an infinitely long cylinder. At the same 
time the longitudinal response shows a dielectric behavior usual for disordered wires. In 
fact, the shorter the localization length in the longitudinal direction is the larger is the 
current in the transversal one. In the limit of high frequencies or short inelastic times the 
transport in the model under consideration is classical in all directions and is described by 
the conventional Ohm law. 

We are grateful to V. Falko for critical reading of the manuscript. 
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